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The Bogoliubov theory is extended to a Bose-Einstein 
condensation with internal degrees of freedom, realized re- 
cently in 23 Na gases where several hyperfine states are si- 
multaneously cooled optically. Starting with a Hamiltonian 
constructed from general gauge and spin rotation symmetry 
principles fundamental equations for condensate are derived. 
The ground state where time reversal symmetry is broken in 
some case and low-lying collective modes, e.g. spin and den- 
sity wave modes, are discussed. Novel vortex as a topological 
defect can be created experimentally. 
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Much attention has been focused on Bose-Einstein con- 
densation (BEC) since its realization in alkali atomic 
gases such as 87 Rb §, 23 Na § and 7 Li @. Furious 
experimental |I| and theoretical [^| investigations have 
been done to reveal various novel phenomena of BEC 
which have not been attainable in another BEC system, 
that, is, superfluid 4 He. One of the most notable differ- 
ences in these two systems lies in their mutual interac- 
tion either very weak in the former or very strong in the 
latter. This enables us to construct a microscopic the- 
ory for BEC in the present system from first principles. 
In fact a mean-field theory or Hartree-Fock-Bogoliubov 
framework due to Bogoliubov |j] , Gross , Pitaevskii ||] 
and others has been quite successful in explaining funda- 
mental physics of BEC realized in magnetically trapped 
atomic gases such as condensate fraction, transition tem- 
perature, or low lying collective modes. We only feed the 
s-wave scattering length a into these microscopic theo- 
ries as a material parameter, which itself is known rather 
accurately j§||]. 

Most experiments for BEC so far were performed by 
the magnetic traps which necessarily select one of sev- 
eral possible atomic ground states, or the so-called weak 
field seeking state, such as F z =-1 among F z — ±1,0 of 
F = 1 in the 23 Na case, where the spin degrees of freedom 
are "frozen" . Recently Stamper-Kurn et al. have suc- 
ceeded in cooling 23 Na atoms purely optically in an opti- 
cal dipole trap and achieved BEC, where the three sub- 
states F z = ±1,0 are simultaneously "Bose-condensed" . 
They demonstrate it by a "Stern-Gerlach" experiment: 
The F = 1 magnetic three sublevels are all visible as a 
split image by passing it through a field gradient. That is 
just remarkable. This opens an interesting possibility to 
explore BEC with internal degrees of freedom where not 
only the gauge symmetry U(l) but also "spin" symmetry 



50(3) for the F — 1 case are involved, a situation similar 
to the superfluid 3 He problem where the orbital (I — 1) 
and spin (s = 1) degrees of freedom give rise to a 3 x 3 x 2 
dimensional manifold of complex order parameter space 

s 

The purposes of this paper are (1) to establish funda- 
mental equations for describing such a situation, leading 
to an extended Gross-Pitaevskii equation, (2) to exam- 
ine low-lying collective modes in order to extract basic 
physical properties of the ground state and finally (3) 
to point out several interesting possible topological ex- 
citations or textural spatial structures under special cir- 
cumstances, e.g. when releasing BEC from an optically 
plugged quadrupole magnetic trap, a novel vortex may 
be created experimentally. So far no one has succeeded 
in observing a vortex even for one-component BEC al- 
though there have been several theoretical discussions 
We will try to demonstrate that the present vec- 
torial BEC system provides us a new horizon and fertile 
field to further explore fll2"f . 

In order to make our argument clear-cut, we consider 
a case described by a constituent atom with the hyper- 
fine state |F| — 1 (F z — ±1,0) (an extension of the 
theory to other cases is easy and trivial) where the or- 
der parameter or the Bose condensate is characterized 
by three components: \I/o, ^-l corresponding to each 
sublevel. Instead of choosing the basis set |1 >, |0 >, 
| — 1 > where the quantization axis is chosen along the 
z-axis, we introduce the basis set \x >,\y >,\z > de- 
fined by Si\i >= (i = x,y,z) where Si is the i-th 
component of the spin operator (\S\ = 1). The order pa- 
rameter is now expressed by a three-dimensional vector 

= (^x, ^y, ^z) on this basis set which behaves as a 
vector under spin space rotation. 

The Hamiltonian density can be constructed which 
should be invariant under spin space rotation and gauge 
transformation (Ti = 1): 
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in terms of the creation and annihilation operators ip^, ip a 
with the a-th species Bose particle (a = x,y,z). The 
Larmor frequency in a vectorial notation ujt, — 7^B (7^: 
the gyromagnetic ratio and B: the external magnetic 
field). The others are the standard notation (Repeated 
index is summed) . The two kinds of the interaction terms 
gi and 32 are characteristic of the multi-component order 
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parameters which represent the spin degrees of freedom 
of the condensate. These are physically originated from 
the interactions: 
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where n — ip^ipa and Sj = ip^Siipa with Si being a 
3x3 spin matrix. It is easy to confirm the relations: 
.<?i — 9n+9s and g 2 = —,9s- Thus the g 2 term in (1) comes 
from the F-spin flipping process. It should be noted here 
that the precise values of g\ and g 2 are neither known for 
each atomic species of interest (Na, Rb and Li) and nor 
determined for each hyperfine state (|F| = 1,2, •••) to 
our knowledge. However, according to a recent spectro- 
scopic experiment by Abraham et al. on 7 Li (Boson) 
which determines the singlet (a s ) and triplet (at) scat- 
tering lengths (a s = 33 ± 2 and a t — — 27. 6 ± 0.5 in units 
of the Bohr length), it is conceivable in general that the 
spin flipping interaction g 2 is comparable to g\ and could 
be either positive or negative. For the present BEC to 
be stable it is necessary that g\ + g2 > as seen shortly. 

Following the standard procedure H] which leads 
to the Gross-Pitaevskii equation and starting with the 
Hamiltonian (1), we can derive the coupled equations for 
the condensate wave functions 
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where in the equation of motion for ip a the ground state 
average < ip a > is replaced by \& Q and the approximation 
such as < ip'tpi/) >^ vj/*ij/\]/ is introduced. This is a basic 
set of equations of a non-linear coupled Shrodinger type. 

When the system is uniform and infinitely large, the 
energy E of the system is given by 



E=- 9l \^ a \ 



:32|* a * a | 2 + ie aPl uj Ll ^* ^ f} . (4) 



The minimization of E is achieved by 



for n g 2 > -uj l 
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with fi = nogi — ujl- The order parameter is intrinsically 
complex and breaks time reversal symmetry, no — ^* ■ $ 
is the density of the condensate. The magnetic moment 
M oc i^P x vjj* whose direction is parallel to the z-axis 
is spontaneously induced in this ground state. This is 
understandable because 172 > or g s < implies the 
ferromagnetic interaction for the spin exchange channel. 
For the opposite case, the minimization is achieved by 
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with fi — n (gi + g 2 ). When B = or (lol = 0), the 
ground state = y/n (l,0,0) is characterized by a 

real order parameter, thus it has no spontaneous moment 
because time reversal symmetry is not broken. 

Since the original continuous gauge and spin space ro- 
tation symmetries are now broken, the Goldstone modes 
must exist. Note that the complex order parameter man- 
ifold is six-dimensional corresponding to the number of 
the degrees of freedom in the present system. In order 
to derive the collective modes, we consider the equa- 
tions of motion for small deviation from the ground state: 
^ a = ^ + S^ a , which are valid up to linear in 5^ a . 
For rio52 > —ll>l, the equations of motion for 8^ x and 
8^ y are coupled, that is, 

d 1 
igl 69 * = ( w i + £ k + My + 92)}S^ X + -nogiS^* 

+i{co L + M~y + 9a)}8*v + ~nogiS% (8) 
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where s k = j— and S^a = S 1 i' a (k, t). These give rise to a 
4x4 eigenvalue matrix while that of 5^ z is detached from 
those. After a little manipulation, we find from Eqs.(8) 
and (9) that the density wave mode and the longitudinal 
spin wave mode are identical, whose dispersion relation 
is given by 



oj = \Js k (e k + 2n .9i) = ck 



with the velocity 



and another mode, 

to = e k + 2(uj l + n g 2 ). 
The equation of motion for 8^ z : 
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where 



leads to the transverse spin wave mode: 

LU = U) L +£ k . (14) 

This massive mode (14) is quadratic in k, corresponding 
to the ferromagnetic-like spin excitation. 

For 77.052 < —tOL similar coupled equations of motion 
give the following collective modes: 
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The former is collective mode coupled with the number 
density and the longitudinal spin density in the presence 
of the external field. The latter is the transverse spin 
wave. When B = 0, these are decoupled to yield the 
density wave mode: 



£? + 2e k n a (g 1 + g 2 ) = ck 



with the zero-sound velocity 

c 2 = "o(ffl + .92) 

m 

and two spin wave modes: 



e\ + 2e k n \g 2 \ = ck 



with the spin wave velocity 
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The former (17) corresponds to the usual Bogoliubov 
zero-sound phonon mode while the latter (19) is the spin 
wave. Both are linear in k in the long wave length limit. 
It is seen from (18) that for the system to be stable 
9i + 92 > is necessary. The order of magnitudes of the 
zero-sound velocity and spin wave one are same. Since 
the former is measured by Andrews et al. Jl5| on 23 Na 
atom gas of a single component BEC, which is an order 
of a few cm/sec and analyzed by the Bogoliubov frame- 
work |l6|-[T8]]. The known scattering length (a=2.75nm) 
turns out to accurately explain the condensate density 
dependence of the zero sound velocity, leading mutual 
reliability between theory and experiment. Here by the 
same token the measurement of the spin wave velocity 
may yield a direct estimate for the other missing impor- 
tant parameter g 2l which determines the physics of the 
system in a fundamental way. 

So far we only consider the spatially uniform case. We 
can expect rich textural spatial structures or topological 
defect structures in this multi-component BEC system 
which is analogous to superfluid 3 He JTcfl . In order to 
help establishing this analogy which may facilitate our 
understanding, it is useful to adopt the following repre- 
sentation: Let us introduce $ = (m + in) with m and 
n real unit vectors and define 1 = m x n. The three real 
vectors (m, n, 1) form a triad and 1 points to the mag- 
netic moment M. Starting with Eq.(l) we can derive 
alternative equations of motion in terms of this triad: 
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where hq = i(m 2 



jg 2 (m 2 - n 2 ) 2 - mV 2 m - nV 2 n 
1 = uj l I x z - Vif s (25) 
n 2 ), the number current 



f p = m • Vjii - n • V,m 
and the spin current 

j* = (Vjin) x m + (Vjn) x n. 
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It should be noted that Eq. (25) does not explicitly in- 
clude the interaction constants g\ and g 2 . When the sys- 
tem is uniform, the I vector precesses with the Larmor 
frequency about the equilibrium position. 





FIG. 1. Possible two dimensional types of the I vector 
texture, (a) radial disgyration, (b) circular disgyration and 
(c) cross disgyration 

The /-vector field could form various types of the spin 
texture jl9), depending on the situation as in superfluid 
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3 He such as the radial, circular and cross disgyrations 
as shown in Fig.l. The last one is particularly inter- 
esting because in this configuration the direction of the 
Z-vector field coincides with that of the quadrupole mag- 
netic field produced by an optically plugged quadrupole 
coil. When releasing a BEC system confined magneti- 
cally by the quadrupole trap, this cross disgyration tex- 
ture may appear as a natural equilibrium configuration. 
Since the order parameter of this texture is described 
by /(r)(sin<^x + cos^y + iz) in terms of cylindrical co- 
ordinate (r, <f>, z), the x- and y-components exhibit the 
circulation ±1, implying that this is a kind of the vortex 
structure. 

It should be noted that the earth field (~0.1G) which 
amounts to the Zeeman energy ~ 10/iK is larger than the 
BEC transition temperature Tg ~ meaning that 

it cannot be neglected. The order of magnitude of the 
permissible residual field after shielding the earth field 
is estimated by equaling the characteristic condensation 
energy griQ ~ f^B ( r o average interparticle distance ~ 
lOOOA and the typical scattering length a ~ lOOA) to the 
Zeeman energy 7 M -B CJ resulting in B c ~ 10 _3 G. This size 
of magnetic shielding is quite feasible. 

The best way to detect the spin excitations is of 
course nuclear magnetic resonance experiment, but in 
view of the small number of the condensed atoms trapped 
(~ 10 6 ) it may not be possible to perform such an exper- 
iment. An optical method utilizing the Faraday rotation 
effect of the polarization change due to the magnetiza- 
tion in a medium may be an alternative. We also notice 
that except for the case of <?2 < under no applied field, 
all the spin wave excitations are coupled with the den- 
sity wave excitations, which can be observed as a phase 
contrast image. 

In conclusion, we have derived the coupled non-linear 
Gross-Pitaevskii equation (3) for BEC with spin degrees 
of freedom, demonstrating that the vectorial BEC is 
much richer internal structures than the spin frozen one- 
component scalar BEC system. In the ground states both 
for <72 > and gn < not only gauge symmetry but 
also spin rotation symmetry are broken. The low-lying 
collective modes, including the gapless Goldstone modes 
are analytically obtained. In order to help facilitating 
discussions on spatio-temporal behavior of the multi- 
component BEC basic equations of motion (21)~(27) for 
a triad (1, m, n) are constructed. The detailed study of 
it belongs to a future work. We propose an experimental 
setup to create a novel vortex by releasing BEC from a 
quadrupole trap. 

We thank D. Stamper-Kurn for useful communication 
on their ongoing experiment and T. Yabuzaki and Y. 
Takahashi for useful discussions. 
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